Thermoelectric properties of Ba-Cu-Si clathrates 
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Thermoelectric properties of the type-I clathrates Ba8Cu3;Si46-i (3.6 < x < 7, x — nominal Cu 
content) are investigated both experimentally and theoretically. The polycrystalline samples are 
prepared either by melting, ball milling and hot pressing or by melt spinning, hand milling and hot 
pressing techniques. Temperature-dependent electrical resistivity, p{T), and the Seebeck coefficient, 
S{T), measurements reveal metal-like behavior for all samples. For a; = 5 and 6, density functional 
theory calculations are performed for deriving the enthalpy of formation and the electronic structure 
which is exploited for the calculation of Seebeck coefficients and conductivity within Boltzmann's 
transport theory. For simulating the properties of doped clathrates the rigid band model is applied. 
On the basis of the density functional theory results the experimentally observed compositional 
dependence of p{T) and S{T) of the whole sample series is analyzed. The highest dimensionless 
thermoelectric figure of merit ZT of 0.28 is reached for a melt-spun sample at 600° C. The relatively 
low ZT values in this system are attributed to the too high charge carrier concentrations. 



I. INTRODUCTION 

Intermetallic clathrates have been investigated ex- 
tensively in recent years as promising thermoelectric 
materials^. They crystallize in several structure types^, 
of which the type I clathrates AgE46 (A = alkali, al- 
kaline earth, or rare earth element, E = group IV ele- 
ment) have hitherto shown the best thermoelectric per- 
formance. The unique features of the clathrate crystal 
structure are generally regarded as ideal realization of 
the phonon glass - electron crystal (PGEC) concept^. 
Clathrates are viewed as Zintl compounds'^'^: the cova- 
lently bonded framework atoms (E) use the valence elec- 
trons of the guest atoms (A) to satisfy the octet rule and 
to acquire closed-shell electronic configurations. Sim- 
ple electron counting based on the Zintl concept^ fre- 
quently works welP'^'^''' for understanding the thermo- 
electric properties of clathrates. Previously, most in- 
vestigations focused on clathrates with heavy framework 
atoms (e.g., Ge- or Sn-based clathrates) which typically 
are partially substituted by other heavy elements (e.g., 
Ga, In, or transition metal elements). Comparably little 
work has been done on Si-based clathrates®"^^, despite 
the technological interest in this low-cost element. Early 
interest in Si-based clathrates was triggered by the dis- 
covery of superconductivity in Na2Ba6Si46 with a critical 
temperature Tc = 4 K (Ref. 16). Later on, investigations 
on Ba8TM3;Si46-x (TM: transition metal element, TM 
= Cu, Ag, Au)^''"^-^ demonstrated the stabilization of 
the crystal structure by TM substitution, allowing for a 
simple synthesis in an arc furnace. However, this led to 
a suppression of which, from band structure calcu- 
lations, could be attributed to a change of the density 
of states near the Fermi leveP^"^'*. Interest in the ther- 
moelectric properties of TM-containing clathrates also 
arose from band structure calculations which indicated 
that Ba8{Cu,Ag,Au}6Si4o clathrates are degenerate p- 



type semiconductors with large Seebeck coefficients at 
room temperature^'^'^^'^^. This motivated us to inves- 
tigate the thermoelectric properties of Ba8Cu2;Si4g_2; in 
detail. 

Recently, we investigated the phase equilibrium of 
clathrates in the Ba-Cu-Si system at 800° C and its 
crystal structure^^. On the other hand, first-principles 
methods have made considerable successes in studies of 
structural and electronic properties of materials. More- 
over, the Boltzmann transport theory has been exten- 
sively applied to investigate electronic transport prop- 
erties such as electrical conductivity and Seebeck coeffi- 
cient. A method that hybridize the first-principle meth- 
ods and Boltzmann theory would no doubt benefit under- 
standing of thermoelectric properties of materials, which 
has proved to be quite successful in our recent work^*. 
In this paper, we present a comprehensive study on the 
thermoelectric properties of Ba8Cu2;Si46-2: (3.6 <x<7, 
X = nominal Cu content) in combination with density 
functional theory (DFT) calculations and the semi-classic 
Boltzmann transport equations. 



II. EXPERIMENTAL 

A. Sample preparation and characterization 

A series of clathrates, with the nominal compositions 
Ba8Cu3;Si46-a; [x = 3.6, 3.8, 4.2, 4.4, 4.6, 5, 6, 7, denoted 
by HP01-IIP08), was prepared in a high-frequency induc- 
tion furnace from high-purity elements (more than 99.9 
mass%), as described previously^''. The samples were 
subsequently annealed (800°C, 15 days), ball milled and 
hot pressed (T = 800°C, P = 56 MPa, time = 2 hours). 
A sample with a; = 6 was chosen for melt spinning'^^''^^ 
with wheel speeds of 1500 r/min and 2500 r/min (sam- 
ples MS1500 and MS2500). 
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X-ray powder diffraction (XPD) data for all samples 
were collected using a Siemens D5000 diffractometer (Cu 
Kai,2, 10° <2e < 110°) and processed with the Rietveld 
method using the program FULLPROF^^. The crystal- 
lite size of hot pressed samples was evaluated from the X- 
ray diffraction patterns using the MDI jade 5.0 software 
(Materials Data Inc., Liverpool, CA). Annealed (300°C, 
1 day) pure Si powder was used for the evaluation of the 
peak broadening due to the X-ray instrument. The bulk 
density of the hot pressed samples was measured by the 
Archimedes method. 

The compositions of the hot pressed samples were de- 
termined by energy dispersive X-ray spectroscopy (EDX) 
in a scanning electron microscope (SEM) operated at 20 
kV (Philips XL30 ESEM; JEOL JSM-T330A, probe size: 
1/um). 



B. Physical properties 

The electrical resistivity and the Seebeck coefficient 
were measured with a ZEM-3 (ULVAC-Riko, Japan). 
The thermal conductivity was calculated from the ther- 
mal diffusivity Dt, measured by a laser flash method with 
a Flashhne-3000 (ANTER, USA), the specific heat Cp 
and the bulk density D using the relation k = DtCpD. 
Hall effect measurements were performed with a standard 
ac technique in a physical property measurement system 
(PPMS, Quantum Design) in the temperature range from 
2 to 300 K in magnetic fields up to 9 T. 



C. Computational details 

Density functional theory (DFT) calculations were 
carried out by the Vienna ah initio Simulation Pack- 
age (VASP)^^'^^. The exchange-correlation functional 
is approximated within the generalised gradient approx- 
imation using the parametrization of Pcrdcw, Burke, 
and Ernzerhof (PBE). The electron-ion interaction is 
treated within the framework of Blochls projector aug- 
mented wave mcthod^'''^^. The valence state configu- 
rations for the construction of the pseudopotentials in- 
cluded 5s^5p^6s^ states for Ba, 3p^3rf^°3s^ states for Cu, 
and 3,s^3p^ states for Si. For all calculations a 5 x 5 x 5 k 
point grid according to Monkhorst and Pack'^^ was used 
to sample the Brillouin zone. 

Transport properties for the electrons were modeled 
within the semi-classical Boltzmann transport theory in 
which the Seebeck coefficient tensor (5) and electrical 
resistivity (p) are defined as 

p = a-i a{T,EF)=e^Ko (1) 



Here Kn is related to electronic structures by 

^" = ^jE^(i'k)v(i,k)(£(i,k) -i;^)" 

dE 

The vectors v = dE{i,'k) /dk represent the velocities of 
the electrons as defined by the derivative of the band en- 
ergy, r, / and Ep are the relaxation time, the Fermi 
function, and the Fermi energy, respectively. Within 
Boltzmanns transport theory effects of electron-electron 
and ck;ctron-phonon scattering arc merged into the re- 
laxation time r. Because a first-principles calculation of 
r is not feasible for systems with large numbers of atoms 
per unit cell, r was considered as an empirical parameter 
by fitting to one selected experimental value at a given 
temperature, as described below. 

It should be noted that the cubic symmetry of Ba8Si4o 
is maintained for Ba8Cu6Si4o but broken when model- 
ing Ba8Cu5Si4i. Therefore, the Seebeck tensor has now 
unequal components. Because the experimental sam- 
ples were produced at higher temperatures the measured 
structures are always cubic with some random site occu- 
pations. In order to enable the comparison with exper- 
iment the symmetry average of the calculated Seebeck 
and conductivity tensor is made, resulting in only one 
component for each physical property. Concerning the 
transport properties of Ba8Cu5Si4i two sets of calcula- 
tions have been made, namely for a geometrically fully 
relaxed structure and for a structure denoted as frozen- 
Ba8Cu5Si4i, for which the high-symmetry structure of 
Ba8CueSi4o was chosen and no relaxation was allowed. 
The actual calculations of the transport properties were 
made by an adapted version of the package BoltzTrap^^. 
For this purpose the Kohn-Sham energy eigenvalues were 
generated on a very dense 25 x 25 x 25 k-point grid. 

III. RESULTS AND DISCUSSION 

A. Sample purity in the hot pressed samples and 
structural details of clathrates 

The phase constitutions of all hot pressed samples 
agree well with the isothermal section of the system Ba- 
Cu-Si at 800°C (the homogeneity range of the clathrate 
at 800°C: 3.4(1) < xcu < 4.8(1), Ref. 27), except for a 
very small amount of Si (and/or of holes) that is uni- 
formly dispersed in the main clathrate phase (even for 
the samples with a: > 5). The compositions from EDX 
measurements are close to the nominal compositions for 
the samples HP01-HP05 which correspond to the single 
phase range of clathrate in the phase diagram. For the 
samples HP06-HP08, a composition of Ba8.iCu4.9Si4i.o 
is determined. Table I summarizes the nominal Cu con- 
tent, the Cu contents determined from EDX and from 
the XRD refinements, as well as the crystallite sizes and 
densities of all specimens. 
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TABLE I. Hot pressed Ba8Cua;Si46-i samples with their nom- 
inal Cu content (Xnom), the Cu content determined from EDX 
measurements (xedx) and from the XRD refinements (xref), 
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Sc (nm) 


HPOl 


3.6 


3.5 


3.61 


92.7 


165(15) 


HP02 


3.8 


3.7 


3.85 


92.8 


125(10) 


HP03 


4.2 


4.2 


4.29 


92.0 


120(5) 


HP04 
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4.72 
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4.9 
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FIG. 1. (Color online) The Rietveld refinement of the powder 
X-ray diffraction data from HP04 shows the main clathrate 
phase with a small amount of Si. The microstructure of this 
sample (insert) shows the main grey phase (clathrate phase) 
with black spots (Si and/or holes). 



Rietveld refinements of the X-ray diffraction data were 
performed with an initial model derived from the single 
crystal X-ray diffraction data of Ba8Cu5Si4i (Ref. 27). 
Here, Ba atoms fully occupy the 2a and 6c sites in the 
crystal structure of the type-I clathrate, Cu atoms share 
the 6d site with Si, and the remaining sites are occupied 
by Si atoms. For this model the Rietveld refinements of 
the data for all samples have high reliabilities and yield 
reasonable structure parameters (see Table II). Unlike 
for the case of several Ge-based TM clathrates^^'^^"^^ , 
here no vacancy could be identified at the Qd site within 
experimental resolution. Figure 1 shows the Rietveld re- 
finement for the sample IIP04 as an example. 
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FIG. 2. (Color online) Temperature dependent electri- 
cal resistivity, p(T), and Seebeck coefficient, S(T), for 
Ba8Cui:Si46-a; (3.6 < a; < 7) plotted as a function of temper- 
ature. The straight lines in S{T) represent linear fits through 
the origin. 



B. Thermoelectric properties 

1. Hot pressed Bus CuxSi46-x samples 

The temperature dependent electrical resistivity, p{T), 
and Seebeck coefficient, S{T), of the hot pressed sam- 
ples show a systematic evolution with the Cu content 
(Fig. 2). All samples show metal- like behavior. The neg- 
ative S{T) values indicate that electrons dominate the 
transport properties. In the Cu-rich samples (IIP06- 
HP08), extrema are observed in p(T) and S{T). 

The overall variation of p{T) and S{T) with Cu content 
can be qualitatively understood within a phenomenolog- 
ical picture of the electronic band structure of TM con- 
taining clathrates^^. In all samples, the Fermi level lies 
within the conduction band. This explains the metal-like 
behavior of p{T) and the negative S{T) values. With in- 
creasing X, the Fermi level moves towards the edge of the 
conduction band. Thereby, the charge carrier concentra- 
tion n decreases and, consequently, the absolute values 
of p{T) and S{T) increase. 

According to band structure calculations for a fully 
ordered structure, Ba8Cu6Si4o should be a p-type semi- 
conductor (see Sect. IIIC and Refs. 25 and 26). The fact 
that only n-type behavior is found in our sample series 
is attributed to the formation of competing phases when 
Xnom exceeds 4.9. This inhibits a further energetical low- 
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TABLE II. Structure data for Ba8CuxSi46-a; (3.6 < x < 7). The data have been standardized by the program Structure Tidy' 



Name 




HPOl 


HP02 


HP03 


HP04 






3.61 


3.85 


4.29 


4.38 


a (nm), Guinicr 




1.03272(3) 


1.03273(5) 


1.03278(3) 


1.03277(6) 


rf = y \Fo - Fj/y Fo 




0.105 


0.122 


0.082 


0.075 


Ri = y\i„- u/y 




0.085 


0.081 


0.073 


0.060 


Bal, in 2a (0,0,0) Beq (B^so 


) 10^ (nm^) 


0.41(9) 


0.62(9) 


0.47(9) 


0.41(7) 


Ba2, in 6c (j,0, \) Beq (Bisc 


,) 10^(nm2) 


1.71(9) 


1.72(9) 


1.63(8) 


1.57(5) 


Ml, in 6d (|,|,0), Occ. 




3.61(7)Cu+2.39Si 


3.85(6)Cu+2.15Si 


4.29(6)Cu+1.71Si 


4.38(5)Cu+1.62Si 


Beq (Biso) 10^ (nm^) 




0.32(8) 


0.43(9) 


0.42(9) 


0.53(9) 


Sil in 16j (a;, x,x), x 




0.1855(3) 


0.1854(3) 


0.1852(3) 


0.1851(2) 


Beq (Bi.o) 10^ (nm^) 




0.54(7) 


0.85(9) 


0.59(9) 


0.66(9) 


Si2 in 24fc {0,y,z), y,z 




0.1198(5), 0.3077(5) 0.1188(4), 0.3091(5) 0.1189(4), 0.3098(5) 0.1192(3), 0.3094(3) 


Beq (Buo) 10^ (nm^) 




0.43(9) 


0.25(7) 


0.43(9) 


0.27(9) 



TABLE III. continued 



Name 




HP05 


HP06 


HPor 


HP08 


Xref 




4.55 


4.72 


4.88 


4.98 


a (nm), Guinier 




1.03280(3) 


1.03286(4) 


1.03285(5) 


1.03287(4) 


rf = j:\Fo -Fe\/j:Fo 




0.107 


0.088 


0.083 


0.088 


Rj^y\io-wyio 




0.086 


0.071 


0.071 


0.065 


Bal, in 2a (0,0,0) Beq (B^c 


,) 10" (nm") 


0.49(7) 


0.41(9) 


0.38(9) 


0.40(9) 


Ba2, in 6c (i,0, i) Beq {Bi,o) W'^{nm'^) 


1.63(6) 


1.55(9) 


1.52(7) 


1.63(8) 


Ml, in 6d (|, i,0), Occ. 




4.55(5)Cu+1.45Si 


4.72(7)Cu+1.28Si 


4.88(6)Cu+1.12Si 


4.98(5)Cu+1.02Si 


Beq (Bi.o) 10^ (nm") 




0.65(9) 


0.68(8) 


0.43(9) 


0.86(9) 


Sil in 16i (a;, x, x), x 




0.1842(2) 


0.1849(3) 


0.1848(3) 


0.1849(3) 


Beq {B,,o) 10'' (nm") 




0.68(9) 


0.62(9) 


0.65(9) 


0.68(9) 


Si2 in 24fc (0,y,z), y,z 




0.1203(3), 0.3096(2) 0.1190(5), 0.3093(5) 0.1200(4), 0.3095(4) 0.1195(3), 0.3092(3) 


Beq (Biso) 10" (nm") 




0.40(9) 


0.52(9) 


0.28(7) 


0.39(9) 



ering of the Fermi level down to the valence band. The 
energy gap is predicted to decrease with increasing Cu 
content (sec Sect. Ill C and Rcfs. 25 and 26). This is in 
agreement with the observed maxima in p{T) and S{T) 
for the Cu-rich samples HP06-08, which signal the onset 
of intrinsic conduction above about 500°C. 

To a first approximation, for temperatures above the 
Debey temperature 9d {Od ~ 150°C for Ba8Cu5Si4i, 
Ref. 47), and ignoring the phonon-drag contribution to 
the total Seebeck coefficient^^, S{T) should vary linearly 
with temperature as^^ 



(4) 



where n is the charge carrier concentration, m,; is the 
free-electron mass, e is its charge and the other symbols 
have their usual meaning. Thus, one should be able to 
estimate n from the slopes of linear fits to S'(T). These 
fits are shown as straight lines in Fig. 2. For high Cu 
contents, these fits describe the data quite well in rela- 
tively broad temperature ranges. For low Cu contents. 



the agreement is poorer. As will be shown below, this can 
be attributed to a temperature dependent charge carrier 
concentration. 

For the samples HPOl, HP02, HP04 and HP06 we have 
also determined the charge carrier concentrations more 
directly, by measurements of the Hall coefficient Rh{T) 
below 300 K. The charge carrier concentrations calcu- 
lated using a simple one-band model nn = ^/{gRh), 
are plotted against temperature in Fig. 3. For HP06, 
nniT) is almost temperature independent. For the Cu- 
poor samples, on the other hand, nniT) shows a sizable 
temperature dependence. 

Finally, we can derive the charge carrier concentration 
expected within a simple Zintl electron counting scheme. 
For Ba8CuxSi4i_a;, the charge carrier concentration per 
unit cell, nz, varies with x as nz{x) = 16 — 3a; , if the 
oxidation states +2, -3, and arc assumed for Ba, Cu 
and Si^** , respectively. A comparison of the charge carrier 
concentration determined via all above methods is shown 
in Fig. 4. 

The temperature dependent thermal conductivity. 
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FIG. 3. (Color online) Temperature dependent charge carrier 
concentration n//(T) of HPOl, HP02, HP04 and HP06 derived 
from Hall effect measurements in a one-band model. The lines 
are guides to the eye. 
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FIG. 4. (Color online) Charge carrier concentrations derived 
from different methods (see text) vs xedx- 



k{T), varies systematically with x, as shown in Fig. 5(a). 
In Fig. 5(b) and (c), we plot the temperature dependent 
electronic component Ke{T) and the lattice component 
K L (T) , respectively. Ke is estimated via the Wiedemann- 
Franz law Ke = LoT/p, where Lq = 2.44 x 10"^ (V/K)2, 
from the electrical resistivity p. kl is the difference 
between k and Ke- The dominating contribution to k 
changes from Ke to k^ with increasing Cu content, with 
almost half to half contributions for the samples HP05 
to HP08. For kl{T) three different characteristics can 
be distinguished beyon the experimental uncertainties, 
as indicated by the gray shaddows shown in Fig. 5(c). 

The Cu-poor samples have the strongest temperature 
dependence and the lowest kl values. With increas- 
ing Xedx and thus decreasing charge carrier concen- 
tration (at elevated temperatures), this dependence be- 
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FIG. 5. (Color online) Temperature dependent total ther- 
mal conductivity k{T) in (a), electronic contribution Ke{T) 
in (b) and lattice contribution kl{T) in (c) for Ba8Cua;Si46-a; 
(3.6 < X < 7). The gray shaddows in (c) indicated different 
characteristics in kl{T) (see text). 



comes more gradual and kl rises. This trend cannot 
result from point scattering or boundary scattering since 
the crystallite size and bulk density do not vary systemat- 
ically with Xedx (Table I). Instead suggest that, besides 
phonon-phonon scattering that usually is the dominat- 
ing scattering at high temperature^^, here also electron- 
phonon scattering contributes sizably. An alternative ex- 
planation for the experimentally observed trend could be 
that the presence of more atoms Cu disturbs the resonant 
scattering of the guest atoms in the large cage. 

The dimensionless thermoelectric figure of merit ZT = 
S^T/pK is plotted against temperature in Fig. 6(a). ZT 
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FIG. 6. (Color online) (a) Figure of merit ZT for 
Ba8Cui,Si46-a; (3.6 < x < 7) versus temperature T. (b) ZT 
vs. XEDx at T = 450°C. 



depends sensitively on xedx, especially for the Cu-rich 
side (Fig. 6). Surprisingly, the samples HP06-HP08 with 
highest X e dx and a small amount of foreign phases show 
the highest ZT values in the investigated temperature 
range. The highest ZT of 0.25 is reached for sample 
HP07 at about 450°C. Between 350 to 550°C, ZT re- 
mains above 0.23. 



2. Melt-spun BasCu^Siio 

In order to further optimize ZT of the clathrates in 
the Ba-Cu-Si system, we melt-spun an as cast sample 
of Ba8Cu6Si4o with two different Cu wheel speeds, 1500 
and 2500 r/min. The resulting ribbons were hand milled 
rather than ball milled to retain the microstructure of 
the ribbons to the largest possible extent, and then hot 
pressed. These samples are denoted by MS1500 and 
MS2500. Figure 7 shows p{T), S{T), k{T) and ZT{T) 
of MS 1500 and MS2500, together with the correspond- 
ing data for the conventionally prepared sample IIP07 of 
the same nominal composition. While the sample IIP07 
reaches the largest absolute S values in the measured 
temperature range, |>5'(r)| of the melt-spun samples is 
still increasing the temperature. The sample MS2500 
has the lowest p in the measured temperature but also 
here a maximum in p(T) appears to be reached at lower 
temperatures for HP07 than for the melt-spun samples. 
MS1500 shows the lowest k{T). Above 480°C, MS1500 
has the highest ZT value. 

As seen from Fig. 7(b), S{T) is very similar for MS1500 
and MS2500. This indicates that also the charge car- 
rier concentrations are comparable: the values 1.88 x 
10^7 m-3 for MS2500 and 1.86 x 10^^ m'^ for MS1500 
are derived from linear fits to S{T) using Eqn. 4 as done 
above. Hall effect measurements also yield similar nn 
within the above used one-band model: 2.18 x 10^'' m"'^ 
for MS2500 and 2.11 x 10^^ m'^ for MS1500. The more 
pronounced differences between MS2500 and MS 1500 in 




7{K) 



FIG. 7. (Color online) Temperature dependent electrical re- 
sistivity, p{T), in (a), Seebeck coefficient, S{T), in (b), ther- 
mal conductivity, k(T), in (c) and figure of merit, ZT{T), in 
(d) of Ba8Cu6Si4o prepared via different preparation routes 
(see text). 



p{T) and k must therefore be attributed to the different 
charge carrier mobilities ph{T) — Rh{T)/p{T), which 
are 2.4 and 3.5 cmVVs at 300 K for MS1500 and MS2500, 
respectively. We attribute that to a smaller average grain 
size in MS1500 than in MS2500. 



C. DFT results 

1. Structure and energetics 

The lattice parameters for the cubic type I clathrate 
structure are derived from fully relaxed VASP calcu- 
lations, from which also the enthalpies of formation 
AH in terms of differences of total energies at zero 
pressure are derived. We compare our calculations on 
Ba8Cu2;Si46-a; with a recent study on the type I clathrate 




FIG. 8. DFT derived formation energies AH and lattice pa- 
rameters of Ba8Aua;Si46-a: (balls) and Ba8Cua:Si46-i (trian- 
gles) as a function of doping. Also shown are the results for 
the unfilled cage trcuture of Si46. 

Ba8Au2;Si46-2;. Figure 8 shows formation enthalpies and 
lattice parameters for four types of systems , namely Si46, 
Ba8Si46, Ba8Cu2;Si46-K and Ba8Aua;Si46-2: with a; = 1 to 
6. From the measured X-ray data it was concluded that 
the Cu atoms are randomly distributed on Si Qd sites due 
to the preparation at elevated temperatures. However, 
for the DFT calculations the doping atoms have to be 
placed on specific sites within the unit cell which, in gen- 
eral, causes changes of lattice symmetry and cell shape. 
Therefore, to be consistent with experiment a cubic unit 
cell was enforced for all calculations. More details are 
given in Refs 28. Comparing to the two experimental 
data points for the lattice parameter in Fig. 8, one finds 
that the DFT derived values are larger by about 1%, but 
show the same trend of increasing with increasing x. The 
(rather small) difference between first principles theory 
and experiment is attributed to the approximation for the 
exchange correlation functional. For the Au-doped com- 
pounds the enthalpy of formation decreases (i.e. bonding 
is enhanced) with increasing Au content, which is very 
similar to the results for Ba8Aga;Ge46-a;^*. The trend 
for the Cu-doped compounds is quite different, because 
ATJ remains rather constant with increasing x. A more 
careful examination reveals that AiJ slightly increases 
for a; > 5 implying that a doping level x exceeding 5 
is unfavorable. This is consistent with the experimental 
fact that no compound can be stabilized for such large 
dopings. 



2. Electronic properties 

The electronic densities of states for Ba8Cua;Si46-a; are 
shown in Fig. 9. Significant changes induced by doping 
Cu are observed. On filling Ba into the voids of the Si46 
cage the gap is maintained but shrinks in size due to hy- 
bridization between Ba and Si states. The gap decreases 
significantly by adding Cu, which implies hybridization 



FIG. 9. Density of states (DOS) of Ba8Cu^Si46-x (a; = 0, 5 
and 6). Fermi energy at iJ = eV. 



between dopants and the Si framework. For Ba8Cu5Si4i 
there are two peaks in the DOS around the Fermi energy, 
which is located at a satellite of the higher energy peak. 
Doping of Cu up to 6 atoms lowers the higher energy peak 
in Ba8Cu6Si4o in magnitude and moves the Fermi energy 
to the center of the peak at lower energy. The small dif- 
ference in the shape of the DOS peaks between these two 
compounds may be due to their different crystal sym- 
metries. However, there is also a significant similarity 
between Ba8Cu5Si4i and Ba8Cu6Si4o in that the peaks 
around the Fermi energy remain almost unchanged and 
also the pseudogap occurs for both of them. The similar- 
ity of the DOS features justifies the rigid band approxi- 
mation when varying the atomic composition in terms of 
changing the number of valence electrons by shifting the 
Fermi energy. For the validity of the rigid band model 
the change of the number of valence electrons must be 
sufficiently small. 

The location of the Fermi level plays an impor- 
tant role in understanding transport properties as dis- 
cussed later. In our recent work within the rigid band 
approximation^*, based on the semiconducting parent 
compounds of Ge46 a simple electron counting rule was 
proposed: Xgap — {n is the valency of the dopant) is 
replaced the Fermi energy in the gap for Ba8M^Ge46_a; 
. This counting rule can also be applied to Ba8M2;Si46-a; 
because of the similarity of Si and Ge. For M — Cu (i.e. 
the valency n equals 1) Xgap = 16/3 is derived, which 
corresponds very well to experimental concentrations at 
which the Seebeck coefficient undergoes a drastic change 
in size. 



3. Transport properties 

For small variations of the dopant concentration (i.e. 
variation An of the number of valence electrons) the rigid 
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FIG. 10. Calculated Seebeck coefficient at 300 K for 
Ba8Cu5Si4i and Ba8Cu6Si4o as a function of doping An. A 
negative value corresponds to hole doping. 



band approximation was employed according to 

/oo 
g{E)]{,T,p)dE (5) 
-co 

for calculating the chemical potential /i (i.e. the Fermi 
energy) for = iVo + Art electrons. The number A^o 
represents the number of valence electrons of the chosen 
reference compound with its DOS, giE). In the Ba-Cu- 
Si clathrates, Ba, Cu and Si contribute 2, 1 and 4 valence 
electrons to the system, respectively. For comparing the 
calculations with measurements on the most interesting 
experimental sample with composition Ba8Cu4.9Si4i.i, 
An = 0.3 if Ba8Cu5Si4i is chosen as the reference and 
An = 3.3 if Ba8Cu6Si4o is the reference. 

Figure 10 depicts the calculated Seebeck coefficients 
at 300 K vs. chemical doping for relaxed- and frozen- 
Ba8Cu5Si4i and Ba8Cu6Si4o, which change sign as the 
Fermi level crosses the gap. There is a significant differ- 
ence in the amplitudes of the Seebeck coefficient between 
the relaxed and frozen structures of Ba8Cu5Si4i reveal- 
ing the influence of structural relaxation. The Seebeck 
coefficient of frozen-Ba8Cu5Si4i resembles quite well that 
of Ba8Cu6Si4o, indicating that the rigid band approxima- 
tion should be useful when varying x between 5 and 6. 
Both curves agree very well with each other except for 
the absolute values nearby the gap. However, the See- 
beck coefficients are significantly different for the fully 
relaxed structures. This is different from the compari- 
son of Seebeck coefficients of Ba8 Ag5Si4i and Ba8 Ag6Si4o 
with their fully relaxed structures, because these practi- 
cally coincide.^* 

Figures 11 and 12 depict the calculated temperature- 
dependent Seebeck coefficients 5'(T) for Ba8Cua;Si46-x, 
which are obviously very sensitive to doping. For com- 
parison, S for Ba8Cu5.4Si4o.6 was also calculated by tak- 
ing Ba8Cu6Si4o as the reference. Although ^(r) for 
Ba8Cu5.4Si4o.6 as derived from the different references 
are qualitatively similar, there are significant differences 



FIG. 11. DFT derived Seebeck coefficients for BagCua;Si46-i 
(x = 5.0 — 5.9) as a function of temperature as calcu- 
lated within the rigid band approximation taking frozen- 
Ba8Cu5Si4i as the reference. For BasCus.4,Si4,os the Seebeck 
coefficient with BagCu6Si4o as the reference is also shown for 
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FIG. 12. DFT derived Seebeck coefficients for BagCua;Si46-a; 
{x = 4.8 — 5.1) as a function of temperature as calcu- 
lated within the rigid band approximation for which frozen- 
Bag Cu5Si4i is the reference. 



concerning the magnitude and the position of the maxi- 
mum. This suggests that care has to be taken for select- 
ing the reference compound. Apparently, S(T) at low 
temperatures changes sign when the Cu concentration 
crosses the critical value of Xgap — 5.33. The change of 
sign for Ba8Cua;Si46-2; can be understood in terms of the 
Mott's formula, which is a simplification valid at very 
low temperatures. There, the sign and magnitude of the 
Seebeck coefficient are determined by the DOS and its 
derivative with respect to energy at the Fermi energy 



dgjE) 



\ep , which results in an opposite sign for 



the Seebeck coefficient and the slope of the DOS at Ep. 
According to Mott a small g{Ep) together with a large 

^^dE ^ |gf gives rise to a large Seebeck coefficient, and 
therefore the Fermi energy should be as close as possible 
to the gap. This consideration is, however, only valid at 
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FIG. 13. DFT derived electrical resistivities of Ba8Cua:Si46-a; 
within the rigid band approximation with frozen-Ba8Cu5Si4i 
as the reference. Crosses denote the positions of the maxima. 



particular position of the Fermi level in combination with 
the temperature dependence of the energy derivative 
of the Fermi function, '^^^^ ■ More specifically, if the 
Fermi energy is above but very close to the gap, the 
states below the gap will contribute to the electronic 
transport only for sufficiently large temperatures when 
^^^^ is sufficiently broad. As the Fermi energy moves 
closer or away from the gap, concomitantly lower or 
higher temperatures are required for thermally exiting 
states below the gap. From the discussion above, one 
derives that the Fermi energies for all the compounds 
are above the gap and are approaching it as x increases 
from 4.9 to 5.1. 

From the discussion above we conclude, that the 
thermoelectric properties and the conductivity of 
Ba8Cu2:Si46-a; alloys can be reasonably described and un- 
derstood by means of a DFT approach. 



sufficiently low temperatures. As temperature increases, 
the Fermi energy moves towards or away from the gap, 
thus changing S'(T). 

From Fig. 9 one derives a positive derivative for 
Ba8Cu5Si4i thus producing a negative Seebeck coeffi- 
cient. For compositions close to Ba8Cu5Si4i, for instance 
Ba8Cu4. 98141.1, the valence electron number is 1.3 elec- 
tron larger than the critical value Xgap = 5.33, which 
means that the Fermi energy is above the gap, higher 
than the Fermi energy of Ba8Cu5Si4i. Further decreas- 
ing the Cu concentration, i.e., increasing the number 
of valence electrons, shifts the Fermi energy away from 
the gap, which significantly reduces the absolute value 
of S{T) although the change in the number of valence 
electrons is rather small. 

Summarizing, the DFT results for S{T) resemble the 
experimental trend rather well. It should be noted that 
rather small uncertainties in the stoichiometry, which are 
unavoidable in the preparation of the samples, result in 
significant variations of the number of valence electrons, 
which concomitantly leads to substantial variations of the 
Seebeck coefficients, as demonstrated by Fig. 12. 

In Fig. 13 the resistivity p{T) is shown for 
Ba8CuxSi46-a; with x = 4.90 — 5.10, as derived for 
frozen-Ba8Cu5Si4i as the reference. The relaxation time 
r = 2.32 X 10"^'' s was chosen to fit the calculated 
p{T) of Ba8Cu5.oSi4i.o to the experimental maximum of 
p{T) of Ba8Cu4.9Si4i.i. Inspecting Fig. 13 one realizes 
that the calculations yield a trend similar to the ex- 
perimental data, because the characteristic peak shifts 
toward higher temperatures with decreasing x. The 
small difference in x between theory and experiment, 
when requiring good agreement, is attributed to the 
experimental uncertainties of sample preparation. The 
maximum in the electrical resistivity is closely related 
to the position of the Fermi energy. In recent work on 
Ba8Aga;Ge46-3;, it was found^* that the peak is due to a 



IV. CONCLUSION 

In summary, thermoelectric properties of the ternary 
clathrate Ba8Cu2:Si46-a; with the nominal Cu contents 
X = 3.6,3.8,4.2,4.4,4.6,5,6,7 have been investigated. 
All thermoelectric quantities vary systematically with 
the actual Cu content in the clathrate phase. The 
highest figure of merit ZT = 0.23 is achieved for the 
sample with the nominal composition cc = 6 at about 
450°C. Almost the same ZT is retained over a rela- 
tively broad temperature range from 350 to 550°C. For 
the optimal Cu content, samples were fabricated also 
by melt spinning. One of these samples reaches an 
even larger ZT value of 0.28 at 600° C. Extensive den- 
sity functional theory calculations have been performed. 
For deriving enthalpies of formation the respective total 
energies for Ba8Si46, ^&^GvlxS\aq-x and Ba8Au2;Si46-x 
were computed. Theroelectric properties were treated 
within Boltzmann's theory and the properties in terms 
of Seebeck coefficients and resistivities were calculated 
by exploiting the electronic structure of Ba8Cu6Si4o and 
Ba8Cu5Si4i. Overall, the agreement between the density 
functional data and experiment is reasonable, but some 
discrepancies remain to be solved. 
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